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Abstract We investigate coefficients in the Israel-Stewart's causal hydrodynamics and discuss the way to calculate 
them microscopic theory. Based on the hadro-molecular simulation based on an event generator URASiMA, we 
evaluate the coefficients for a hot and dense hadronic fluid. 
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Introduction 



A hydrodynamical model is one of the most estab- 
lished models for the multiple-production in high 
energy reactions. Especially for RHIC, a hydrody- 
namical model is believed to be the promising model 
for the V2- Up to now, most hydrodynamical mod- 
els for RHIC adopt a perfect fluid model, i.e., Euler 
equation neglecting viscosity and heat conductivity 
for simplicity. 

It is well know that the relativistic extension of 
the Navier-Stokes equation which contains trans- 
port coefficients is not unique. At least three types 
of deferent definition for the local four velocity J7'' 
are known. In the textbook by Landau and Lif- 
shitz, is defined as an eigenvector of the energy- 
momentum tensor and it is named e-frame. A 
mass current J'' is used to define by Eckart and 
Namiki-Iso^(this is named m- frame). Weinberg 
used a charged current instead of a mass current^. 
All these three types are the covariant extension 
of the non-relativistic Navier-Stokes equation which 
belongs to the diffusion type parabolic equation. 
The rank of the spatial derivative and the one of 
time derivative are different. This difference be- 
tween space and time contradicts relativity and can 
lead acausal results. 

Israel and Stewart proposed a relativistic causal 
hydrodynamic equation'^ which is hyperbolic type 
and contains several additional coefficients other 
than ordinary transport coefficients. In this paper 
we will investigate the physical meaning of these 



new coefficients and evaluate them of a hadronic 
gas based on a microscopic theory. 



2 Causal Hydrodynamics 

The relativistic hydrodynamical equation proposed 
by Israel and Stewart is as follows; 
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The most important feature of Israel-Stewart's 
hydrodynamics is the existence of the relaxation 
terms, (3i, which makes the equation hyperbola. As 
the same order term, there also appear terms with 
ai which stand for the influence of the derivative of 
the different currents. 



3 Linear Responce Theory 

A hydrodynamical model is a macroscopic phe- 
nomenological model of which coefficients should 
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be derived by a microscopic statistical physics. Ac- 
cording to the hnear response theory, transport co- 
efficients are evaluated from the correlations of the 
corresponding currents. Thcrmo-dynamical forces 
and currents are identified through the formula of 
entropy production. 

Israel and Stewart gives entropy production as^l, 
-dad^a - didaa^;)} 
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This formula is diagonal with currents and there ex- 
ists no combination which contradicts with Currie's 
theorem. Therefore, even there exist a term, terms 
added by Israel and Stewart are able to be rewrite in 
the higher derivative of thermo-dynamical quanti- 
ties. Under the basic assumption of linear response 
theory that the disturbance of macroscopic current 
is small, we can safely adopt usual Kubo formula 
for the transport coefficients, such as viscosity and 
heat conductivity [5]. 



4 Evaluation of /3 

As we have discussed in the previous section, we 
can calculate viscosity and heat conductivity by 
taking current-current correlation as usual. If we 
use Israel-Stewart's causal hydrodynamics as a phe- 
nomenological model, we need additional coeffi- 
cients a and f3 also. Neglecting a for simplicity, 
let us focus our discussion on (3 through out this 
paper. 

If the P belongs to a kind of material constant, it 
should only depend on the thcrmo-dynamical quan- 
tities such as temperature and density. It should 
not depend on the boundary condition of the macro- 
scopic current. Hence, we may assume particular 
situation without generality. 

Suppose the situation that the temperature gra- 
dient only exists. Then thermal current q^ in eq. 
(1) is given as. 



If the heat current relaxes in as, 

q\t)^q^O)e-^. (4) 
Then, f3i are related to the heat conductivity as, 

K01 = T^. 

According to the linear response theory, heat con- 
ductivity is calculated through the correlation of the 
heat currents, 

K=i Jd'x'j' dt'e-<'-''\q,ix,t),q,ix',t')). 

(5) 

Therefore, if the heat current relaxes in the relax- 
ation time Tk, heat conductivity is written as, 

K^^«q\t)q\t)»^ 

where << >> is average and integraion in 
cq.(5). 

Hence, transport coefficient and corresponding (3 
are closely related through the expectation value 
of the square of the current, relaxation time and 
temperature. 



5 Relaxation of hadron gas 



• E = 0.370 GeV/fm 
O E = 0.463 GeV/fm' 
■ E = 0.556 GeV/fm' 
□ E = 0.648 GeV/fm'' 

♦ E = 0.741 GeV/fm' 
V = 1 0° f « E = 0.833 Gevnrn 

= 0.231 [fm"'] 
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Figure 1 : Relaxation in visous-shear tensor correla- 
tion 

In a previous paper we have reported the calcu- 
lation of transport coefficients of hadron gas based 



on a hadro- molecular simulation . We used event 
generator URASiMA as a time evolution genera- 
tor. In our simulation, we first put only baryons 
in the box with periodic boundary condition and 
turn switch on the simulator. Collisions between 
particles take place and mesons and resonances are 
produced. At later time, the system approaches to 
some kind of stationary state where all kinds of par- 
ticle possess a common slope parameter, which we 
may call "temperature"^. The configurations of 
the stationary state we use as statistical ensembles. 

Figure 1 displays correlation of stress-shear ten- 
sors as a function of time interval of the currents, 
which correspond to {Txy{x,t),Txy{x' ,t')) in the 
formula of the linear response theory. Figures 2 and 
3 display the shear viscosity and heat conductivity 
of hadron gas obtained by our simulation. Accord- 
ing to our simulation, both transport coefficients 
depend on baryon number density very weakly and 

ternperature dependences are as strong as almost 
2^5 [6] 

On the other hand, the changes of the relaxation 
times of the currents are very mild as functions of 
temperature (figs. 4 and 5). Most part of the tem- 
perature dependences of the transport coefficients 
come from the change of the intensity (expectation 
value of the current square (figs. 6 and 7)). 
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Figure 4; relaxation of visous-shear tensor 
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Figure 5: relaxation of heat conductivety 
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Figure 2: Shear visosity 



V = 10 fm 
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Figure 3: Heat concuctivity 



6 Concluding remarks 

In this paper we have investigated the coefficients of 
the causal hydrodyanamics and evaluated them by 
using a hadro-molecular simulation. According to 
our simulation based on URASiMA, the relaxation 
time of stress-shear tensor and heat current of the 
hadron gas are about 2 fm. The changes with tem- 
perature and baryon number density are very small. 
This results show rather clear contrast against the 
relaxation time in diffusion process which changes 
clearly with temperature ^ ^ . 

The importance of quantitatively investigating 
the relaxation in order to justify the macroscopic 
model has been pointed out by Iso, Mori and 
NamikiM. 
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Figure 6: The square expectation of heat current 
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Figure 7: The square expectation of visous-shear 
tensor 



tremely helpful. 
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